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THIRD-ORDER CONTRIBUTIONS TO ELECTRICAL 


CONDUCTION IN PLASMAS 

By Willard E. Meador 
Langley Research Center 

SUMMARY 

The electrical conductivity is calculated for a highly nonequilibrium plasma corre- 
sponding to the application of an electric field varying linearly with time. Two methods 
are employed: (1) the direct approach consisting of an exact solution to a modification of 

the third-order kinetic equation of Chapman and Enskog, and (2) Everett's technique for 
closing out the macroscopic equations of change with the Grad 13-moment velocity dis- 
tribution function. Comparisons between the two sets of results indicate substantial 
third-order differences for most interparticle interaction potentials of practical interest; 
hence, the Grad 13-moment approximation does not appear to describe adequately certain 
higher order contributions to plasma transport coefficients. 

INTRODUCTION 

Highly nonequilibrium plasmas corresponding to large electron currents are of 
considerable interest in such potential applications as high current pinches, arcs, and 
discharges and in the description of phenomena occurring near electrodes and in low- 
pressure discharges. Departures in the formulation of problems of this type from those 
satisfying the ordinary linear flux theory are usually reflected in the dependence of 
transport coefficients upon the electron diffusion velocity. 

Examples of the more rigorous treatments of high current plasmas include the use 
of Grad's 13 -moment velocity distribution functions by Everett (ref. 1) and Yen (ref. 2) in 
closing out the macroscopic equations of change for the number densities, diffusion 
velocities, heat fluxes, and pressure tensors. Such procedures yield relations between 
the 13 moments and the applied force fields and also between the moments themselves, 
the accuracies of which depend in large measure upon the ability of the Grad approxima- 
tion to predict the various collisional transfer terms adequately. Previous investiga- 
tions of this ability have been performed by Meador (ref. 3) for first-order relations 
between the heat flux, the entropy density, the entropy density production rate, and the 
electron diffusion velocity, and for second-order contributions to the electron pressure 


tensor. Similar tests of Everett’s and Yen’s higher order corrections to the electrical 
conductivity have not been reported. 

The purpose of the present research is to derive the high velocity electrical con- 
ductivity of a simple plasma from an exact solution to a modification of the third-order 
kinetic equation of Chapman and Enskog (ref. 4). Comparisons with similar calculations 
employing the aforementioned Everett procedure yield the following results: (1) signifi- 

cant differences occur in the contributions to the current density which are proportional 
to the cube of the electron diffusion velocity, and (2) the reliability of the Grad 
13-moment approximation is strongly dependent upon the interparticle interaction poten- 
tials in that it is especially poor for very soft and very hard molecules. The second 
result is understandable in light of the fact that the Grad approximation is an expansion 
(ref. 5, p. 22) in the eigenfunctions of Boltzmann's binary elastic collision operator for 
Maxwellian molecules; consequently, a rapid convergence for more general force laws is 
not guaranteed. 

The particular plasmas chosen for this study are those for which the following con- 
ditions are applicable: the electron collisions obey the collision model recently devel- 
oped by Meador (ref. 6), the heavy particles are infinitely massive and at rest relative to 
the laboratory, the applied and induced magnetic fields are zero (the latter in strict vio- 
lation of Maxwell's equations, but consistent with a nonrelativistic treatment), all macro- 
scopic quantities are spatially homogeneous, and the applied electric field is proportional 
to the time. 


SYMBOLS 

a electric field parameter defined in equation (3) 

A ratio of collision integrals defined by equation (21) 

b impact parameter 

c e electron particle velocity relative to laboratory frame of reference 

c h heavy -particle particle velocity 

e magnitude of electron charge 

E applied electric field 

f e electron distribution function 
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Maxwellian contribution to electron distribution function 
heavy-particle distribution function 
,gg trial functions appearing in equation (18) 

index numbers 
electron current density 
Boltzmann's constant 
unit vector in z-direction 
electron particle mass 
electron number density 
number density of heavy particles 
electron partial pressure 

traceless electron pressure tensor relative to laboratory frame of reference 

traceless electron pressure tensor relative to electron frame of reference 

integral defined by equation (10) 

entropy density 

equilibrium entropy density 

collisional production rate of entropy density 

time 

electron temperature corresponding to zero electron diffusion velocity 
electron temperature relative to laboratory frame of reference 
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electron temperature relative to electron frame of reference 

reduced electron particle velocity relative to electron frame of reference 

electron diffusion velocity 

Cartesian coordinates; also used as subscripts to indicate vector and tensor 
components; in addition, x is used as an integration variable 

electron body force per unit mass 

third-order electrical conductivity parameter 

reduced electron diffusion velocity 

reduced heat flux relative to electron frame of reference 

reduced electron particle velocity relative to laboratory frame of reference 

azimuthal angle for collisions 

conductivity parameter defined by equation (55) 

interparticle interaction parameter 

electrical conductivity for zero electron diffusion velocity 
effective electrical conductivity through third order 
collision time for electron diffusion 
collision time for entropy production 

first-, second-, and third-order electron perturbation functions, respectively 
scattering or deflection angle 
collisional time derivative of f e 
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Primed quantities in collision integrals refer to conditions after a collision; 
unprimed quantities, before a collision. When vector symbols appear without an arrow, 
the magnitude of the vector is denoted. The symbol ( ^ indicates an average over 
velocity space. 


KINETIC THEORY 


The derivations of high-velocity electrical conductivities and other plasma trans- 
port coefficients are usually based in some manner upon the electronic Boltzmann equa- 
tion (ref. 4) 



( 1 ) 


where y is the reduced electron particle velocity defined by 



(2) 


and T® is the electron temperature at zero electron diffusion velocity. 

If the time -dependent electric field and electron velocity distribution function are 
expressed as 

E = kat (3) 

and 



respectively, the following simplification of equation (1) is obtained from the application 
of the special conditions outlined in the Introduction and the resulting stipulation from the 
electron equation of continuity that n e is constant: 


5 



( 5 ) 


More specifically, the successive simplifications of the collision integral in equa- 
tion (5) correspond to the assumptions that the heavy particles are fixed scattering cen- 
ters for the electrons and that electron-electron collisions can be explicitly neglected at 
this stage. The latter assumption is part of Meador's collision model (ref. 6), the 
remainder of which outlines the method whereby the single electron — heavy -particle 
interaction potential can be generalized semiempirically to include the effects of multiple 
heavy species as well as electron-electron encounters. These more detailed aspects of 
the collision model will appear subsequently in the form of an effective interaction 
parameter which mathematically (but not physically) assumes the role of the exponent 
in an inverse-power electron — heavy-particle interaction potential. No restrictions are 
yet placed upon the characteristics of the unknown functions <p 2 > and ^3. 

The method adopted in the present research for the direct solution of equation (5) 
involves the following three statements: (1) the velocity distribution function is analytic 
in the small electric field parameter a, (2) the unknown function <j>^ } which is also small 
and satisfies <pi > contains only the ith power of a, and (3) the index number i of 

4>i designates the order of solution. Hence, the retention of functions through (£3 
implies a third-order solution in the sense that contributions to f e of all a 3 terms are 
included. Of particular significance in this procedure is the fact that the corresponding 
first-order form of equation (5), which is written as 



differs from that of Chapman and Enskog (ref. 4) in the explicit appearance of <p^ on the 
left-hand side. 
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Equation (6) is further simplified by the assumption that <p ^ is y z multiplied by 
a function of y and possibly the time. The following equation results 


8 ^1 + 2eat/ m e = _ R 13 

8t m e UkT°y z r 04 t o 



when the T°-dependent collision time 


(7) 


and the collision integral (ref. 6) 



bdbde = 


_ m e CT o 


e 2 n e 


R 13 j 

/m e V 

n h R 04 T o( 

V 2kT°y 


r 7? 


( 8 ) 


(9) 


are employed. The integrals are defined by 



As mentioned previously and explained in detail in reference 6, the semiempirical 
interaction parameter £ in equations (7), (9), and (10) appears formally as the exponent 
in an inverse-power electron — heavy-particle interaction potential and can be so chosen 
as to make the present theory reliable for many real plasmas (electron-electron colli- 
sions included). Only in the case of Lorentz plasmas (slight ionization or fully ionized 
gases with large ionic charges), however, can £ be physically identified with such an 
exponent, the value of which may range from unity (Coulomb forces) to infinity (rigid 
spheres). Even if the generalization is not possible, the present analysis would still 
accomplish its primary purpose of evaluating the capability of Grad 13-moment distribu- 
tion functions to describe highly nonequilibrium situations. 


An exact solution of equation (7) is readily found to be 


*1 


2eaR Q4 T Q / me 
m e R 13 \2kT° 


a/2 


R 04 t o £ 

y 

r 13 



(ID 


7 


so that the first-order current density becomes 



if r s is given by 

^ r 04 r 22 ^ 

T c = T, 



(12) 


(13) 


The r s of equation (13) is shown in reference 6 to be never less than r 0 nor 
more than twice t 0 and to satisfy the entropy-production definition 

\ T S s c = s(0) " s ( 14 ) 

Accordingly, t s is called the characteristic time for the collisional production rate of 
entropy and appears as a convenient parameter in many applications; for example, the 
kinetic theory expression of the Hall conductivity for small magnetic fields assumes the 
simple mean-free-path form if r s is substituted for t 0 . The present interest lies, of 
course, in the fact that r s acts as a lower limit for t in the expression for the cur- 
rent density, and thereby signifies the size of the time scale inherent in the Boltzmann 
collision integral. 

As a final comment on the first-order solution, it is immediately obvious from 
equation ( 11 ) that an upper limit on the time (and thus the electric field) must exist if 
( p 1 is to be small compared with unity. This requirement, however, does not create any 
special problems prior to the consideration of $> 3 , and will be discussed again at that 
point. 


SECOND-ORDER SOLUTION 

The substitution of equation (11) and its velocity derivative 
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into equation (5) yields 
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£y 2 - 2^y 2 + | - 4)y| ^ = -n h y^_ ®j j* ^> 2 + ^3 " ^2 “ 0;j) b db de 


(16) 


The second-order form of equation (16) is obtained as follows by deleting all terms 
corresponding to a 2 : 


90 2 e2a 2 R 04 T Q t | “ 3 f 

3t + m e kT°R 13 ? y 


£y 2 - 


4 _ j 

2£y 2 + £ - 4)y 2 - - y^ £y 2 - 2^y 2 + £ - 4)yi 


= - n h rlrs^j j '( 02 -^) bdbde 


,o\ 1//2 


If a second-order function of the type 


( 17 ) 


^2 


gi(r) + g2(y)t + g3(y)^ 


i^y 2 - y 2 ) + g 4 (y)t 2 + g 5 (y)t 3 (18) 


is assumed, the following equation results: 
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by using the collision integral 




b db de = 
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2n h R 04 T ol2kT° / / 




derived in reference 6. 


The quantity A is the ratio of collision integrals 


A = 



- cos 2 y 


b db 



- cos x)b db 


Equation (19) is automatically normalized in the sense that 


j 4 O) 0 2 dc e = 0 

In addition, the temperatures T e and Tg relative to the laboratory and electron 
frames of reference, respectively, are given by 



where is the magnitude of the reduced electron diffusion velocity defined as 



It is further noted that the use of equation (12) in the time derivative of equa- 
tion (23) yields the principle of conservation of energy 




= E • j 


for this problem; thus, all the pertinent auxiliary conditions relative to n e and T e 
satisfied by the $2 of equation (19). 


(20) 


( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 


(26) 

are 
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Because sizeable values of t q J t are not especially important, it is convenient 
now to examine the second-order solution after T Q j t has become small compared with 
unity. The neglect of squares and higher powers of this time ratio in equation (19) gives 


2e 2 a2R 4 r 0 t 3 T" 1 
02 = r r 
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(27) 


(28) 


(29) 


as the components of the traceless electron pressure tensor P e relative to the labora- 
tory frame of reference. 

Since the traceless electron pressure tensor relative to the electron frame of ref- 
erence has a zz-component given by 

4p J3 2 
> _ e 
ezz 3 


o T 

p ezz = n e m e \P-ez 


the following equation can be written: 


P T 

r ezz 


; - v e > 

) ~ n e kT; 

4p e /S 2 

16(| + l)T g 
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15A|t q 


- 1 


from equations (29) and (30) if T s j t is neglected. In particular, 

% zz (Z = 1) = l-42p e ^ 2 

as compared with the Grad value of 1.17p e /3 2 from reference 3. 


(30) 


(31) 


(32) 
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THIRD -ORDER SOLUTION 


The magnitude of T Q /t is next assumed to be the same order as the reduced elec- 
tron diffusion velocity /3 for the purpose of solving the third-order kinetic equations. 
Since t q times the third term on the left-hand side of equation (16) can be written as 


with the aid of equations (12) and (25), and since the next to last term in the <fi 2 of equa- 
tion (19) similarly becomes 


9 4 

4Rq 4/3 t ^ 

9£ r 13 t o 



the combination of these two expressions is proportional to (fit/ r 0 and is therefore of 
order in magnitude and third order in the electric field parameter a. The corre- 
sponding contributions from the remaining terms of <p 2 have orders of magnitude ($, 
;3 4 , and /3 5 and are neglected in the present section. 


It is further evident from the contribution to (p 2 which is retained, and which 
includes the factor t Jt q , that some such establishment of the order of magnitude of 
Tq/ t (and hence the introduction of an upper limit to the elapsed time and the applied 
electric field) is necessary in order for <p 2 to be kept smaller than (pi and thus for 
the present expansion to be convergent. This aspect was previously anticipated. 


The appropriate reductions in equation (19) and its velocity derivative for use in 
equation (16) are 


and 


0 2 - 


2e 2 a 2jD T +3 4 
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9 7z 
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(33) 


(34) 


Thus, the third-order kinetic equation becomes 


12 



(35) 


8 ^3 + 2e 3 a 3 R Q4 T o t 4 /na { 

8t 9m^kTeRi3l 2 \ 2kT ^ 




2£ 2 y 4 - |(3| + 8)y 2 - (| + 2)(| - 4) 


/2kT 


oN 1 / 2 


= ~ n hA 


m c 


J (<#>3 ' 03) 
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Equation (35) can be solved in a similar manner to equation (6) to yield 

4R2 4% ^Et/ me \ 1 / 2 f-4 | 


0 3 =_ ' 


9en e R 2 ? 2 r 0 \2kT° 


13 


2£ 2 y4 - 5(3^ + 8)y 2 - (| + 2)(| - 4) 


y y 


(36) 


through the leading contribution. Accordingly, if the current density is expressed as 


j 




1 + aU) 



(37) 


the electrical conductivity parameter a is derived from the velocity moment of equa- 
tion (36) to be 


4(g + 2) g - 4)r s t 
HU + 8)r 2 


(38) 


A brief analysis of the essential features of this entire development shows the 
leading term of the fifth-order contribution to the current density to be proportional to 
j3 5 (t/r 0 ) , which is third order in magnitude if r Q J t and j3 are regarded as first order 
in magnitude. Equations (37) and (38) thus represent the complete description of j 
through terms which are third order in a and second order in magnitude. They also 
comprise the exact third-order result mentioned in the Introduction as a purpose of the 
present research. 


THE GRAD APPROXIMATION 

The final effort of the present research is the determination of a using the Grad 
13-moment approximation with Meador's collision model (ref. 6) to close out the macro- 
scopic equations of change. Numerical comparisons with equation (38) should yield valu- 
able insight into the applicability of this technique to highly nonequilibrium plasmas 
because a completely common framework (that is, basic assumptions and expansions in 
powers of the electric field parameter a) is provided for both methods. Hence, the 
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preceding exact solutions are unquestionably the correct standards for such comparisons 
and the only critical question is whether the 13 -moment function adequately satisfies the 
pertinent kinetic equations or their moments. 

Although not necessary, a direct confirmation of the statement that the preceding 
exact solutions are the correct ones is obtained from reference 3. Additional moments 
caused the Grad approximation to converge rather rapidly toward similarly derived exact 
solutions in that application. 

Since Everett's velocity distribution function (ref. 1) can be written as 


f e = n e 


m. 


,3/2 


-u 


27rkTi 


1 + * u2 " 2 K ‘ * + ( n e kT eJ p e : 


-lo 


(39) 


when the reduced heat flux and the reduced electron particle velocity relative to the elec- 
tron frame of reference are defined by 


and 



(40) 


(41) 


respectively, the following third-order closed-out expressions are obtained from the 
m e c ez and m e c ! c ez moments of equation (1) and the use of equations (9) and (39) in 
evaluating the collision integrals: 
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e 


/ \ 1/2 

(| - 4)/3’ /2kT'\ 

+ ±- — ® =0 

5£ \m e J 


(42) 


and 



10Rj^tcx o /3 2 E 
9en e Ri3R_i,5(3£ - 2 )t q 


(43) 


As in the preceding exact development which culminated in equations (36) to (38), all but 
the leading ^tfr 0 terms have been deleted from these expressions. 

The simultaneous solution of equations (42) and (43) finally gives 


5U " 4)R 2 ta oi 3 2 E 

v P (third order) = — — 

9en e (3| - 2)RjgR_^ 


(44) 
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so that 


j (third order) = -en e Ve(third order) = - 


5(1 - 4)R^ 4 ta oi 32 
9(3£ - 2)R 13 R_ 1 5 t 0 


E 


(45) 


A direct determination of a from equation (45), however, is not possible unless 
the Grad approximation also yields equation (12) through first order in a. In particular,, 
adjustments may have to be made in order to account for differences between the predic- 
tions of cr 0 by the exact and Grad methods. This calculation starts with the following 
first-order closed out expressions analogous to equations (42) and (43): 


d/3 eat / m e 


m e V 
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and 


d/3j 
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dt m ey2kTg 
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5Ro 4 ! 2 t 0 


5^/3 + (3^ - 4)/3j 


(47) 


J3 - - 


and 


The simultaneous solution of equations (46) and (47) gives 

R 04 R 13( 179 ^ 4 ' 496 ^ 3 + 832 ^ 2 " 512 £ + 256 ) t c 


R 2 4 (l3^ 2 - 164 + 16)a 0 E/ m X 1 / 2 


4en e Ri 3 R_i )5 (3^ - 2)4 \2kT°, 
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5R 0 S - 4 )°o E he 
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x R 04 Ri 3 (234 2 -164 + 16)t £ 


4R 22 R-i 5 (34 - 2)4t 


(48) 


(49) 


Accordingly, equations (45) and (48) combine to yield 


- rL(i 34 2 - 164 + leWr 
j= 7 — ^ l+o<4) 

4R 13 R -1 5 ( 3 ^ - 2)1 L 


R 0 4 R 1 3 ( 4 7 9 4 ^ - 4964 3 + 8324 2 - 5124 + 256)t, 

4R 22 R -1,5 (3 ^ ' 2 )( 13 ^ 2 “ 16 4 + 1 6)4 t j 


and 


a = 


204(4 - 4)t 


9(l34 2 - 164 + 16) T 0 

if proper regard is taken of the order of magnitude consistent with this development. 


(50) 

(51) 
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Percentage errors in the ratio of Grad to exact conductivities are obtained from 

°o(Grad) _ R 04( 13 ^ 2 ~ 16 ^ + 16 ) (c 

Ob 4R 13 R_ 1 j 5 (3? - 2)g 

and range from -4.3 at g = l throughout g = 4 to -4.3 at g = As mentioned in 
reference 3, on the other hand, a somewhat different picture prevails for the heat flux 
because equation (49) must be compared in that case with 



/ \l/2 / 

(g - 4)o 0 E/m e \ / 

2en e^ \2kT°y \ 



(53) 


from the energy moment of equation (11). For example, Grad's result is in error by 
-26 percent at g = 1, if T g /t is neglected when compared with unity. 


COMPARISONS OF RESULTS 


A convenient summary of the preceding research on electrical conduction (through 
second-order terms in magnitude and third order in a) is provided by the expression 

J = crE (54) 

where o is the effective electrical conductivity defined by 

a = a Q (exact or Grad)^. - + a/3 2 ) (55) 

The parameter t) assumes the values 

T/(exaet) = 1 (56) 

from equation (37) and 

„ (Grad) = *0 (57) 

4R 22 R -1,5 (3 ^ " 2) ( 13 ^ 2 " m + 16 ^ 

from equation (50), whereas the formulas for a are given in equations (38) and (51). 

Numerical calculations appropriate to equations (54) and (55) are presented in 
tables I and II for a variety of effective interparticle interaction potentials ranging from 
the fully ionized Lorentz plasma (g = 1) to a gas of rigid spheres (g = °°). Only in the 
neighborhood of Maxwellian molecules (g = 4) is the Grad 13-moment approximation 
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adequate for the determination of a, the error being very significant for both softer and 
harder force laws. The differences however are more tolerable in the case of rj. 


TABLE I.- THE ELECTRICAL CONDUCTIVITY PARAMETER a AS A FUNCTION 
OF THE EFFECTIVE INTERPARTICLE INTERACTION PARAMETER | 

AND THE METHOD OF SOLUTION 


s 

0 tT o jt (Grad) 

ar 0 jt (exact) 

1 

0.513 

0.859 

2 

.247 

.196 

4 

.000 

.000 

OO 

-.171 

-.524 


TABLE H.- THE ELECTRICAL CONDUCTIVITY PARAMETER rj AS A FUNCTION 
OF THE EFFECTIVE INTERPARTICLE INTERACTION PARAMETER £ 

AND THE METHOD OF SOLUTION 


1 

?](Grad) 

rj(exact) 

1 

0.759 

1.000 

2 

1.011 

1.000 

4 

1.000 

1.000 

oo 

.861 

1.000 


CONCLUDING REMARKS 

Calculations through third order in the electric field have indicated that large 
errors can occur when the Grad 13-moment velocity distribution function is used to close 
out the macroscopic equations of change. Except in the neighborhood of Maxwellian force 
laws, the only transport coefficient (among those considered) for which the Grad approxi- 
mation yields accurate results is the first-order electrical conductivity. The difficulties 
begin with the first-order heat flux relative to the electron frame of reference and the 
leading contribution to the traceless electron pressure tensor, each of which is in error 
by 20 to 30 percent for very soft or very hard interaction potentials, and are magnified 
many times in the case of third-order transport phenomena. Although the plasma chosen 
for the present work is a hypothetical one, this trend of the Grad 13-moment approxima- 
tion toward greater (and problem-dependent) discrepancies seems to be well established 
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for highly nonequilibrium systems; consequently, some of the past research on large 
electron diffusion velocities should perhaps be reevaluated. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 19, 1969. 
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